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Abstract — Logistic functions are good models of biological 
population growth. They are also popular in marketing in 
modelling demand-supply curves and in a different context, to 
chart the sales of new products over time. 

Delays being inherent in any biological system, we seek to 
analyse the effect of delays on the growth of populations governed 
by the logistic equation. In this paper, the local stability analysis, 
rate of convergence and local bifurcation analysis of the logistic 
equation with one and two delays is carried out and it can be 
extended to a system with multiple delays. 

Since fluctuating populations are susceptible to extinction 
due to sudden and unforeseen environmental disturbances, a 
knowledge of the conditions in which the population density is 
fluctuating or stable is of great interest in planning and designing 
control as well as management strategies. 

Index Terms — delay logistic equation, stability, rate of conver- 
gence, Hopf bifurcation analysis. 

I. Introduction 

The logistic equation is a simple model of population 
growth in conditions where there are limited resources. It was 
proposed by Verhulst in 1838 to describe the self limiting 
growth of a biological population [1]. The equation has a 
variety of applications. It is used in neural networks to clamp 
signals to within a specified range [2], in economics to 
illustrate the progress of the diffusion of an innovation through 
its life cycle [3], in medicine to model the growth of tumours 
[4] and in linguistics to model language change [5]. 

The growth of a population is conceptualised in Fig. 1. 
A biological population with plenty of food, space to grow 
and no threat from predators tends to grow at a rate that is 
proportional to the population. However, most populations are 
constrained by environmental limitations. Growth is eventually 
limited by a factor, usually one from amongst many essential 
resources. When a population is far from its limits of growth 
or the carrying capacity of the ecosystem, it can grow ex- 
ponentially. The feedback about the availability of resources 
reaches with a delay due to various factors such as generation 
and maturation periods, differential resource consumption with 
respect to age structure, hunger threshold levels, migration 
and diffusion of populations, markedly differing birth rates 
in interaction species and delays in behavioural responses to 
a changing environment (including changes in density of prey 
or predators or competing species) [6-8]. When nearing its 
limits, the population can fluctuate, even chaotically. 

According to the logistic equation, the growth rate of a 
population is directly proportional to the current population 



and the availability of resources in the ecosystem. The logistic 
equation model is as follows, 




where y is the population at that instant, a is called the 
Malthusian parameter which represents the growth rate and 
k is the carrying capacity of the ecosystem. 

Hutchinson incorporated the effect of delays into the lo- 
gistic equation [9]. Delays bring about interesting topological 
changes in the population size like damped oscillations, limit 
cycles and even chaos [10]. The bifurcation analysis of a 
system with a single delay has been performed in [11]. The 
importance of two delays in the logistic equation can be seen 
in [12-15]. In this paper, the effect of such delays will be 
analysed. Analysis will primarily focus on local stability, local 
rate of convergence and local bifurcation phenomena. The 
procedure adopted here to characterise a system with two 
delays can be easily extended to a case of multiple delays. 

The rest of the paper is organized as follows. In section 
2, the model is described. In section 3, It is linearised, local 
stability analysis of the model is carried out and the conditions 
for stability are presented. In section 4, rate of convergence 
is analysed. Bifurcation analysis of the system is carried out 
in section 5. Finally in section 6, the graphical results of the 
analysis are presented. 

II. Model Description 

In this section, we present the model employed with the 
accompanying assumptions for a single delay system and a 
two delay system. 

A. Logistic equation with single delay 

The model we have employed to demonstrate population 
dynamics is as follows: 

• The initial normalised population is chosen to be small 
(typically 0.01) as it cannot be zero. A zero initial 
population signifies a non-existent species. 

• The equation is normalized. 

• The growth rate is a which is finite, positive and time 
independent. 

• The delay t\ is also finite and positive. 

With these assumptions, the single delay logistic equation 
simplifies to 
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Gestation. 
Maturation periods, 
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Figure 1. Schematic of the model. The population growth depends both 
on the current population (reproductive growth) as well as the availability of 
resources. The logistic equation with delay model the abstraction well and is 
widely used. 



^ = ax(t) [1 — bx(t — ti)] . 



(2) 



B. Logistic equation with multiple delays 

The same assumptions as made is the single delay logistic 
equation my be applied here as well. In the case of a system 
with n delays, the equation takes the form 



dx 
~dt 



ax(t) 



1 - YjbiX(t - Ti) 



(3) 



The case of two delays is analysed in this paper. 



III. Local Stability Analysis 

The delay logistic equation is non-linear. In this section, 
we first linearise it and proceed to extract conditions for the 
stability of both the single delay as well as the multiple delay 
logistic equation. We end the section by presenting stability 
charts of the system. 



1 ) Stability analysis around initial value: Substituting x* 
or linearising around the initial value, we have 



dp 
~dt 



ap(t), 



(5) 



the solution of which is easily obtained as 

p(t) = he at , 



which is to say 



x(t) = be at , 



(6) 



around x = 0. This equilibrium point is not stable as the 
characteristic equation of (5) has as its root a which is in the 
right half plane. 

2 ) Stability analysis around saturation value: This value is 
of greater interest to us. Substituting x* = \ or linearising 
around the saturation value, we have 



dp 
dt 



-ap(t - ti). 



To solve this, we take 



p(t) = b 2 e xt . 
Then, the characteristic equation reduces to 



A + ae 



-An 



0. 



(7) 



For the system to be stable, the roots of the above equation 
must lie on the left-half of the A plane. Bifurcation point is 
the value of the parameters for which the roots lie on the 
imaginary axis. 

Substituting A = juj in (7), we get 

juj + ae~ jujTl = 

=>• juj + a [cos(co>ti) — j sin(a;Ti)] =0. (8) 
Equating real part to zero, we get 

a cos(cjti) = 



A. Stability analysis of logistic equation with single delay 

Let the equilibrium point (where the growth rate is zero) be 
denoted by x*. Then x(t) = x(t - ti) = x*. Also, 



^ = = ax*(l 
dt v 



bx* 



(4) 



Solving for this, we get x* = (signifying the initial stages) 
and x* = | (signifying the saturating stages). 

Now, to linearise the equation, substitute x(t) = x* +p(i), 
where p(t) is a small variation in the population. Higher 
powers of p may be neglected. The equation therefore is 



dx 
~dt 



= a [x* + p(t)] [1 - bx* - bp(t - n)] . 



=^un = (2ra+l)-,rc = 0,1,2... 
Equating imaginary part of (8) to zero, we get 

uj — asin(u;Ti) = 0. 
To get minimum order solution, we use n = in (9). 

7T 

2-ar 1= 0. 
For r = 0, the characteristic equation is 



(9) 



A 



-a, 



which implies that A < for a > and the system is stable, 
i.e., for ar\ = 0, the roots of the characteristic equation 
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lie of the left-half plane. Therefore, necessary and sufficient 
condition for stability is, 



7T 

an < -. 



(10) 

It is apparent that if the system has no delays or t\ = 0, the 
system is always stable. 



B. Stability analysis of logistic equation with two delays 

At the fixed points, x(t — r±) = x(t — r 2 ) = x(t) = x*. 
Here x*is the equilibrium point. The delay logistic equation 
then becomes: 



dx 
~dt 



= ax* [1 - (&! + 62)2;*] . 



(11) 



Equating ^| to zero, we get to be or bl ] rb2 • We have seen 
in the previous subsection that the first equilibrium is unstable 
and is not of much interest. We explore the latter value. 

The delay logistic equation is then linearised by expanding 
it by a Taylor series and neglecting all higher order terms. We 
obtain the following linearised equation about the fixed point 



6i+6 2 



, such that x(t) 



y(t). We obtain: 



dy 
dt 



= -ax* [hy(t - ti) + b 2 y(t - r 2 )) . 



(12) 



To solve this, we take 



y(t) 



he xt . 



If bi ^> 62, then the second delay term in (12) can be neglected 
and the analysis reduces to the case of a single delay system. 
Consider the case when b\ = b 2 , 



Ae At = _^ e A(t-r l)+e A(t-r 2 )^ 

The characteristic equation therefore is 

A + V Ari + V Ar2 = 0. 

For the system to be stable, the roots of the above equation 
must lie on the left-half of the A plane. Bifurcation point is 
the value of the parameters for which the roots lie on the 
imaginary axis. Substituting A = juo in (13), we get 



(13) 



juj + ■ 



+ -e~ juJT2 = 



ju + - [cos(cjti) - j sin(u;Ti)] 

+ ^ [cos(cjt 2 ) - j sm(ujr 2 )] = 
Equating imaginary part of (13) to zero, we get 



0. 



cos uo — - — cos uj — - — = 0. 



The second term cannot be as the condition on the imaginary 
component will imply that uo = which is not true for all uo. 



2(±f + 27rCi) 
cu — , 

where C\ G Z. 

Using the value of u obtained in (14), we get 
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fl "^+^)(™(^)-™(^))' 

where (n + r 2 ) [sin(^) - sin(^)| + 0. 
On simplifying this, we get the bifurcation point at 



a = 



(14) 



(15) 



(16) 



(-l+^COs(f^)' 

The sufficient and necessary condition is: 

a(r 1 +r2)cosf ^ (ri ~ T2 j ") < tt. (17) 

Substituting n = r 2 = r, we get the familiar single delay case 
ar < |. These conditions are valid for any positive value of 
a, nand r 2 . We also get the following sufficient condition for 
the specific case of b\ = b 2 from (17), 

+ T 2 ) < 7T. 

C. Sufficient conditions for stability 

In this sub-section, we explore the sufficient conditions for 
stability in both the single as well as multiple delay models. 
Consider the characteristic equation of a system with two 
delays: 

A + ax% ie - Xri + ax*b 2 e~ XT2 = 0. 
This can be re- written as follows [16], [17]: 



*r *r \Ti<2X*bi / \ v 

A + ax*bi + ax*b 2 H ^- (e~ Xri - l) 

An 



^r 2 ax*b 2 ^_ Ar2 _ ^ 

AT2 



We now define 



Hi(X) = A + ax*b\ + ax*b 2l 



and 



uo — - sm(cjri) — - sm{UT 2 ) = 0. 



Equating real part to zero, we get 



a a 

- COS{UJTi) + - COS((jJT 2 ) = 



e~ Xri — 1 \ / e" 

#2 (A) = Ariax*6i ( — — ) +Ar 2 ax*6 2 



Aro 



Ar 2 



Clearly, i^i (A) has no zeros on the right half plane. Consider 
the imaginary axis, 
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ffi(A) |>| A | . 



For any real 0, 



1 - e~ Jb 
I- — — 



< i. 



Hence, on the imaginary axis we have, 

| H 2 {\) \<\ A || riax*6i + r 2 ax*b 2 



If 



then: 



riax*6i +T2ax*62 < 1, 



H 2 {\) \<\ A | . 



(18) 



Hence, by Rouche's Theorem, i^i(A) + H 2 (X) ^ on the 
imaginary axis. There cannot be any zeros in the right half 
plane and the system is stable. Hence a sufficient condition 
for stability is: 



ajbiTi + b 2 r 2 ) 



< 1. 



(19) 



For a system with only one delay, i.e. b 2 = 0, the sufficient 
condition is 



an < 1. (20) 
It is noted that the condition (19) is conservative. 



D. Stability charts 

In Fig. 2, the stability chart of a system with single delay 
is seen. The plot simply corresponds to the bifurcation point 
condition ar = |. Fig. 3 has the stability chart of the two 
delay system. Fig. 4, the stability of the two delay logistic 
equation is shown with respect to the growth rate and the two 
delays. 



E. Nyquist Plots 

In this subsection, we see a few representative Nyquist plots 
of the system when stable or unstable. 

1) Single delay system: In Fig. 5, we observe the Nyquist 
plot of the characteristic equation (7). If there are any en- 
circlements about the origin, the system is unstable as the 
characteristic equation does not have poles. 

2) Two delay system: In Fig. 6, the Nyquist plot for a 
system with two delays is shown. There are no encirclements 
about the origin, hence the system is stable. This matches with 
our calculation as the parameters satisfy the conditions in (17). 
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Figure 2. Stability chart of a single delay system. Neither the feedback 
delay nor the growth rate must be too large for the system to be stable. The 
region that is stable is below the necessary and sufficient condition and Hopf 
bifurcation occurs on this line. 
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Figure 3. Stability chart of a two delay system. Parameter a = 1. The region 
below the curve is stable and corresponds to the necessary and sufficient 
condition for stability. 



IV. Rate of Convergence 

In this section, the analysis of rate of convergence is carried 
out for a system with a single delay. We define the rate of 
convergence (R) as the inverse of the settling time based on 
a tolerance band of ±36.8%. 

Consider the characteristic equation (7). If there are no 
oscillations in the system, the imaginary part of A is zero. 
Let 5ft {A} = a, where 5ft(x) denotes the real part of x. The 
equation now becomes, 
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Figure 4. Stability chart of a two delay system. The system is stable below 1 i I I I r~ 

the surface and enters a Hopf bifurcation at the contour shown. It is not locally 5 10 15 20 

stable for points which lie above the contour. 
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Figure 5. Nyquist plot for an unstable system with a single delay, a = 
2, r = 1 in this diagram and it does not satisfy the sufficient and necessary 
conditions highlighted in (10). Encirclements are seen about the origin. 



g + ae 



= 0. 



On differentiating (21) with respect to a, we obtain, 



1 



0. 



Eliminating a from (21) and (22), we find, 

1 



(21) 



(22) 



(23) 



This point can be considered as the critical damping point. 
For values of ar < ^,the system behaves in an overdamped 
fashion and converges without any overshoot. When \ < 
ar < |, underdamped behaviour is observed, i.e. convergent 
oscillations are present. 



Figure 6. Nyquist plot of a system with two delays, a = b\ = 62 = 1. 
r\ = 0.5 and T2 = 1 is chosen such that conditions in (17) are met. There 
are no encirclements about the origin and the system is stable. 



The characteristic equation (7) can be rewritten as in [18], 

Arie Ari = -an 



A 



W(-an) 



n 



where W(x) is the Lambert W function. The rate of conver- 
gence R is given by, 



R 



3? 



W(-an) 



(24) 



Rate of convergence charts 

In Fig. 7, the rate of convergence of the logistic equation 
with single delay is seen with respect to the time delay. It 
can be seen that the system converges fastest for a slightly 
underdamped system. Fig. 8 shows the variation of rate of 
convergence of solutions to (7) with respect to the growth rate. 
It can be observed that rate of convergence tends to zero as 
growth rate tends to zero. In other words, a population which 
multiplies very slowly requires a very large time to reach its 
carrying capacity. Fig. 9 shows the variation of the rate of 
convergence with both parameters of the single delay system. 

V. Bifurcation Analysis 

In this section, it is shown that the delay logistic equation 
undergoes a Hopf bifurcation at a critical value of the parame- 
ters growth rate (a) or delay (ti). The nature of the bifurcation 
is characterized in the subsequent sub-sections as are the 
properties of the resulting periodic oscillations. Bifurcation 
diagrams are presented at the end of the section. 
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Figure 7. The rate of convergence is plotted with respect to the time delay 
in a system with single delay where the growth rate a = 1. It is seen to 
converge fastest for a slightly underdamped system when riis slightly greater 
than the critical damping point r± = -. 
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Figure 8. The rate of convergence is plotted with respect to the time delay 
in a system with single delay where the delay n = 1. 



A. Existence of Hopf bifurcation 

From (10), at the critical point for the single delay case, 
clti = |. Differentiating characteristic equation (7) with 



respect to the growth rate, we get: 



dX 
da 



-An 



ar\e 

Evaluating at critical point a* = 



-An 



2n ' 



3? 



da )a=a* 



> 0. 



(25) 
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Figure 9. Rate of convergence with respect to growth rate a and delay 
r\ in the single delay system. The red portions represent the fastest rate of 
convergence and blue the slowest. The white portions are for the regions 
where the system is not stable. As can be seen from the chart, fastest rates 
of convergence are obtained for systems with large values of growth rate and 
a small value of the delay. 



Similarly differentiating (7) with respect to the delay , we 
obtain: 



d\_ 



Xae~ Xri 



1 — ar\e 
Evaluating at critical point t{ 



-An * 



2a' 



3? 



dX 1 

^Jn: 



> 0. 



(26) 



The transversality condition of the Hopf spectrum with 
respect to the growth rate and time delay is satisfied in (25) 
and (26) respectively. Thus, logistic equation with single delay 
undergoes a Hopf bifurcation at the critical point given by 

an = f . 

Consider the logistic equation with two delays given by 
(13). Differentiating (13) with respect to the growth rate we 
get, 



dX 
da 



-An 



-AT2 



arie~ Xri + ar2e~ Xr2 



Evaluating at the critical point given by 

a 



(T 1+ T 2 )C0S($^)' 



(27) 



we get, 



3? 



{-} 

[da ) c 



a(ri + r 2 ) + 7T 

-2 + a 2 r 2 + 

1 2 1 ri+r 2 



> 0. 



(28) 



Hence the transversality condition of the Hopf bifurcation 
is satisfied with respect to the growth rate. Thus, the two delay 
logistic equation undergoes Hopf bifurctaion with respect to 
the growth rate at the critical point given by (27). 
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B. Direction and stability of the Hopf bifurcation in the 

logistic equation with growth rate as the parameter u t (t) = A^u t + R^u t . (34) 

The logistic equation with a single delay or with two delays xhe b if urcating solutions of u(t, /i(e)) of (29) have amplitude 
undergoes a Hopf at a critical value of the growth rate as 0(e)? period P(e) and non _ zero Floquet exponent /3(e), where 
shown in the previous sub-section. In this section the direction, ^ p and p have the following expansions: 
stability and period of the bifurcating solutions is analysed. 

The procedure adopted is based on the centre manifold theory ^ = ^ 2e 2 ^ ^ 4e 4 _j_ 0(e 6 ) 

[19] (See also [20]). The logistic equation with two delays as 
described in (3) can be written as, 

P = 4n (1 + T 2 e 2 + T 4 e 4 + 0(e 6 )) , 

ii(t) = ax(t) x / (x(t - ri),x(t - r 2 )) , 
where, /3 = /3 2 e 2 + /3 4 e 4 + 0(e 6 ). 

If \i 2 > 0, the bifurcation is supercritical and if /i 2 < 0, it 

/ (x(t -r 1 ),x(t-r 2 )) = l- (hx(t - n) + 6 2 a;(* - r 2 )) . is subcriticaL If & < 0, iz (*, /x(e)) shows asymptotic orbital 

^29) stability and instability if /3 2 > 0. These coefficients will now 

Without loss of generality assume that r 2 >r 1 . (29) can also be calculated - 

be written as Define the adjoint operator A* as, 

u(t)=L»u t + F(u t ,ri, (30) f-f, se(0 , Tl ] 

= t > 0, /J, E R, where for r 2 > 0, U-r 2 *7 s = 0. 

Here, i(j e C ([0, r 2 ], 5ft)and r] T denotes the transpose of rj. 

ut(0) = u(t + 0),u: [-r 2 ,0]^R,0 G [-r 2 ,0]. For ^ G C ' ([-^2,0], 5R) and ^ G C([0,r 2 ],») define an 

inner product 

Also, L M : C[-t 2 , 0] is 

/•0 /»(9 

^0 = -(oo + (6i0(-n) + M-r 2 )) <V>, 0) = ^(0)0(0) - / / ^ T (f - 0)^(0, 0)0(0^. 

J0=-t 2 J%=0 

and F(ix t , /i) : C[-r 2 , 0] R is (36) 

Let </(0) be the eigenfunction for A corresponding to A(0), 
namely 

F(u u fji) = -(ao+fijx* (biu t (0)u t (-ri) + b 2 u t (0)u t (-r 2 )) . 

Here a = a + /i and a refers to the critical value of the A o q(0) = juj q(6), 

growth rate at the bifurcation. 

By the Riesz representation theorem, there exists a q(@) = e Ju; ° 6 \ 

matrix-valued function with bounded variation components 

r,(0, /i), G [-r 2 , 0], such that * Let 9*00 be an ei § en vector of A o such that 



M= / d77(0,M)0(0), (3D 

«/ — T2 

where G C ([-r 2 , 0], 3?) and, 



and 

<</*,</> = 1, ((7*,g)=0. 



dry(0, /i) = -(a + /x) x [6x5(0 + n) + + r 2 )] . From the above equation D can be determined as shown 

below. 

Here 5(6) is the Dirac-Delta function. 

For cj) G C 1 ( [-72 , 0] , 3?) , define o 



*,q) = D-D [ [ e- juJ ^- e Ur](e)e juj ^d^ 
Je=-T 2 Je=o 



A 6(9)- i ^ 0G[-r 2 ,0) 
" m ~\l°. T2 dv^mO = L^ = 0, (32) 

and ^l = D-D J 6e juj ° e d^ff) 



-T 2 



I e [-r 2 , 0), ^ 1 = £)_ D - X * [ Tie -j»oTi (ao + + T2 e-^or 2(ao + ^^J 

f 9 = 0. " i 

=> L> = 

Now, (30) can be written as 1 - (o + [JineJ^ + 6 2 T 2 e^<^] 
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It can be easily verified that (q*,q) = 0. 
Define, 

z(t) = (q*,u t ) , 
w(t,0)=ut(9)-2to{z(t)q(0)}. 
Then on the centre manifold Co, 



j(t,0)=w(z(t),z(t),0), 



where 



z 2 z 2 
w(z,z,0) = w 2 o(#)y + wn(6)zz + w 02 (6)y + 

and z and z are local coordinates for centre manifold Co in 
the direction of q* and q*. Note that w is also real if u t is 
real, we consider only real solutions. For solutions u t G C 
of (11) at /i = 0, 



£(*) = (q*,Aut + Rut) 

= juj z(t)+q*(0)F (z,z) 

= ju z(t) +g(z,z), (37) 



ut(0) = w(z, z, 0) + zq(6) + zg(<9) 

2 -2 
= W2o(0)y+Wll(0)*^ 

from which i^(0) and u t (—ri) can be obtained. As we only 
need the coefficients of z 2 , zz, z 2 and z 2 z, we get 

ixt(0) = w(z, z, 0) + z + z, 



txt(-n) = wfo ^, -n) + ze~ ju;ori + ze J ' 



^t(0)^t(-n) = w(0)w(-ti) + w(-r 1 )(z + 2) 



+w(0) ( 



- r i _|_ ze^^ 1 ) 



+ zz (e jU0Tl +e-^ 0Tl ) + z~e 



ix t (0)ix t (-r) = 



2? 2 e -j Wo ri + z 2 e juJoTl + ^(e ju;ori + e" ju;ori ) 



+z 2 z(2^n(0)e-^ ori + 



^20 (0) 



-J^ori 



+W 11 (-T) + ) + ... 



here, 



Now we can write, 



z 2 z 2 z 2 z 

g(z,z) = q*(0)F (z,z) = g 20 — +gnZZ+go2 y +#21 ~Y+- 

Now using (34) and (35) we get, 

w = u t — zq — zq (38) 



or, 



#20 



<t(0)F (z,z) 

z z z z 

920 y +5n zz +#02y +521-^- + ••• 



-g*(0)2ax* [he^ 1 + b 2 e juJ0T2 ] 
-max* [6ie- jW0Tl + & 2 e~ jw ° T2 ] (44) 



w 



Aw-2$l{q*(O)F o q(0)}, 



0e[-n,O) 



Aw-2M{q*{0)F q(6)} + F , = 0, 
which can be written as 



w = Aw + H(z, z,0), 



(39) 



gn = -Dax* [6i(e ju;oTl + e - j "° Tl ) + b 2 {e ioj ° T2 + e"^ 07 " 2 )] , 

(45) 

502 = -D2ax*[b 1 e j " 0Tl + b 2 e j "° T2 }, (46) 



where 



521 



2 -2 

H(z,z,0) = H 2O (0)^+H 11 (0)zz + H O2 (0)^ + ... (40) 

Expanding the above series and comparing the co-efficients 
we get, 



-D2ax*^b(w 11 (0)e- juJOTl + ^^l e joj ° Tl 
+wii(-n) + W2 ° ( ~ Tl) ) + b 2 (w u (0)e-^ 



)} 



(2juj o -A)w 2o (0)=H 2O (0), 



-Awu=H n (0), 



-{2jLv o + A)w O2 {0)=H O2 {0). 
From (38), we have 



(41) 



(42) 



(43) 



We have, for e [-n,0), 



H(z,z,0) 



2m*{0)F q{6) 

-gq{0) - M(0) 



z z 

920 ~~~T + 9ll*Z + 902-pr + • 



2; £ 
520TT + flll^ + fl02TT + ■ 



q{0) 
q(0) 
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which yields 



20 



~92oq(0) -go2q(0), 



Hii = -giiq(0)-giiq(O). 
Using the above in (41) to (43), we get 

w 2 o(0) = 2juj w 2 o(0) +g2oq(0) + So2q(0) 



where c\ (0) is the lyapunov coefficient. The asymptotic form 
of the bifurcating periodic solutions is 



u(t, /i(e)) = 2eJ? (q(0)e juJot ) + e 2 J? (E ie 2juJot + £ 2 ) + 0(e 3 ) 

(61) 

for < t < P(e). 

On substituting n = r 2 in (29) the analysis reduces to that 
(47) of the single delay logistic equation. 



wn(0)=g 11 q(0)+g 11 q(0). 
On solving the above differential equations, 



W 2 (0) 



wn(0) 



g2o 



q(0)e j 



gp2 



q(0)e 



(48) 



(49) 



— ^OW" " - — q(0)e- joJot> 
JU jwo 



(50) 



for some £"i , E 2 which will soon be determined. 
Similarly, for 9 = 0, 

ff( 2 ,2,0) = -3?r-F 1 o9(0)+-F 1 o 



ff 2 o(0) 



-520?(0) - 5029(0) - 2a(b ie - ju ^ 

+6 2 e-^ 0T2 ) (51) 

ax* (b 1 w 2 o(-Ti) + b 2 w 20 (-r 2 )) (52) 

+2juj w 20 (0) (53) 



Hn(0) 



-5ii9(0) - 5n9(0) - 2a &i(e 



ax* (6iWn(-ri) + b 2 w 11 (-T2 j) 



JJ0T1 ^ 

(54) 
(55) 



The expression for tU2o(— n)> u> 2 o(-r 2 ), w>2o(0), wii(-ti), 
wn(— T2), and wn(0) can be found from (49) and (50). On 
substituting them in (51) and (52), E\,E 2 can be found. From 
the above analysis the following can be calculated: 



ci(0) = f 520511 - 2|5n| 2 - ^ I502I 2 



521 

2 ' 



M2 



»(ci(0)) 



»(A'(ao))' 
P = 4r(l + e 2 T 2 + 0(e 4 )), 

^ = _ / 9M0|+M2A>q)) ' 
V 

/3 = e 2 /? 2 + C>( e 4 ), /3 2 = 2»( Cl (0)), e = 



(56) 

(57) 
(58) 

(59) 
(60) 



C. Direction and stability of the Hopf bifurcation in the 
logistic equation with delay as the parameter 

In the previous sub- section the direction and stability of 
the bifurcation was considered with the growth rate as the 
bifurcation parameter. It has been shown in Section 5.1 that 
the single delay equation undergoes Hopf bifurcation with an 
increase in the delay also. A similar analysis for this case is 
carried out here. 

Let r be the critical value of the delay at bifurcation. The 
logistic equation can also be written as, 

ii(t) = L^u t + F(u t , /i), (62) 
where, u(t) = xfai), t > 0, /i G R, n = r + /i, and 

u t (0) = u(t + 6)u : [-1,0] -> R,0 e [-1,0]. 
Also, L M : C[-1,0] -> R is 



Lpcj) = -ab(f)(-l). 
F{u u p) : C[-1,0] is 

F(0,/x) = -aty(O)0(-l). 



(63) 



(64) 



By the Riesz representation theorem, there exists a matrix 
function with bounded variation components 77(0, /i),0 G 
[-1,0], such that 



(65) 



where, 



dry(0,/x) = -a6£(<9 + l)d0, 



where <5(6 > )is the Dirac delta function. 

Proceeding further as shown in the previous sub-section we 
get, 



(3tt - 2) 

Hence the Hopf bifurcation is supercritical. 



(66) 
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I 

0.0 0.5 1.0 1.5 2.0 2.5 
Population 

Figure 10. The state space diagram for the two delay system as shown in 
the numerical example with the value of a = ao + 0.05. 




1.15 1.20 1.25 1.30 

Growth Rate (a) 



Figure 11. The bifurcation diagram for the system with two delays as 
described in the numerical example. As can be seen, the results obtained 
through simulation match the analytical predictions for values of the bifurca- 
tion parameter fi = a — ao < 0.1. 

Numerical Example: If we take the case of a system with 
a two delays with b\ = b 2 = 0.5, t\ = 1 and r 2 = 2, 
the system undergoes a Hopf bifurcation at growth parameter 
ao = 1.2092. The value of the lyapunov coefficient ci(0) = 
—0.1691 — 0.2290z, the real part of which is less than zero 
which makes the resulting periodic solutions asymptotically 
orbitally stable. Parameter /i 2 = 0.5175 > 0, implying that 
the Hopf bifurcation is supercritical. The period of these 
oscillations is P = 6. Fig. 10 has the resulting state space 
diagram where the limit cycles can be seen and Fig. 11 
compares the bifurcation diagram obtained through simulation 
to the analytical solution. 

D. Bifurcation diagrams 

In Fig. 12 and 13, the bifurcation diagrams of the single 
delay system are drawn with respect to both parameters 




o 

o 

n 1 1 1 

0.8 0.9 1.0 1.1 1.2 

Growth rate (a) 

Figure 12. The bifurcation diagram of a single delay system with respect to 
the growth rate a. Here the value of the delay t\ = ^ . As can be evidenced 
from the graph, the equilibrium value is constant with respect to a before the 
critical value a = 1. The amplitude increases sharply with changes in a. 




o 
o 



0.8 0.9 1.0 1.1 1.2 

Delay (t<i) 

Figure 13. The bifurcation diagram of a single delay system with respect 
to the time delay r\ while the growth rate a = ^. The equilibrium point 
is not dependent on n when the system is stable and the amplitude of the 
periodic oscillations vary as much with variations in t\ beyond the critical 
value n = 1 as it did with variations in a. 



namely, the growth rate a and delay t\. Both curves look 
similar. In Fig. 14 and 15, the bifurcation diagrams of the two 
delay system are shown with respect to various parameters. 
The only parameter that affects the equilibrium value is the 
scaling factor of the delay terms, b\ and b 2 . The amplitude 
of the periodic oscillations beyond the critical value of the 
parameter varies more with changes in the growth rate a than 
it does with delays t\ and r 2 . 




11 




0.2 0.4 0.6 0.8 1.0 1.2 1.4 



bi 

Figure 14. The bifurcation diagram of a two delay system with variations in 
62 when a = b\ = r\ = 1 and T2 = 10. The equilibrium value depends on 
this parameter. The amplitude of periodic oscillations does not vary as sharply 
with variations in 62 beyond the critical value. 








0.0 0.5 1.0 1.5 2.0 

Delay (x 2 ) 

Figure 15. The bifurcation diagram of the two delay system with parameter 
T2 when a = b\ = 62 = 1 and t± = 2. The variation of amplitude of 
periodic solution is gradual with variations of T2 beyond the critical value. 



VI. Results 

In this section, we show the behaviour of the delay logistic 
system with time under various conditions. In Fig. 16, we 
see the sigmoid curve of the logistic function. As growth rate 
is increased, the function reaches saturation value faster. In 
Fig. 17, the time domain response of a converging system is 
observed. In Fig. 18, the logistic equation with two delays is 
observed in the time domain. Varying the parameters changes 
the behaviour from asymptotically stable to forming limit cy- 



o 




w ' 1 1 1 1 1 ' 

-4 -2 2 4 

Time(t) 

Figure 16. Time domain response of logistic equation with no delay (a = 1). 
Note the sigmoid shape of the curve. The initial growth rate is exponential 
but as the population approaches the maximum carrying capacity, it saturates 
due to diminished resources. The logistic equation was a correction to Robert 
Malthus' proposition that population growth is purely exponential. 




5 10 15 20 25 

Time(t) 

Figure 17. Varying coefficient of the growth rate term a at constant value 
of delay n = 1. This may be viewed as the population charts of different 
species, each with a specific resource consumption rate per capita. The species 
that consume more are the ones that saturate at a lower population. Also, the 
behaviour of a stable and converging system can be seen. 

cles. In the latter case, the sufficient and necessary conditions 
for stability are not met. 

VII. Conclusions 

We have performed local stability analysis of the logistic 
equation with and without multiple time delays. Sufficient 
conditions to aid design were also extracted. The rate of con- 
vergence for the single delay logistic equation was analysed. 
While the ideal system (without delay) is always perfectly 
stable, the actual system that we have considered using delays 
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Figure 18. Varying T2 at r\ = 1, second delay coefficient being twice the 
first, i.e. 62 = 26i. When T2 = 2 or lower, the system is converging to an 
equilibrium point. When T2 = 10, limit cycles are seen. 



undergoes a Hopf bifurcation which is supercritical. The 
nature of the resulting periodic oscillations was analytically 
characterized for the single delay system and a methodology 
to ascertain the same in the two delay case was presented. 
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